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The role of the charm quark in the dynamics underlying the AI = 1 /2 rule for kaon decays can be 
understood by studying the dependence of kaon decay amplitudes on the charm quark mass using 
an effective AS = 1 weak Hamiltonian in which the charm is kept as an active degree of freedom. 
Overlap fermions are employed in order to avoid renormalization problems, as well as to allow 
access to the deep chiral regime. Quenched results in the GIM limit have shown that a significant 
part of the enhancement is purely due to low-energy QCD effects; variance reduction techniques 
based on low-mode averaging were instrumental in determining the relevant weak effective low- 
energy couplings in this case. Moving away from the GIM limit requires the computation of 
diagrams containing closed quark loops. We report on our progress to employ a combination 
of low-mode averaging and stochastic volume sources in order to control these contributions. 
Results showing a significant improvement in the statistical signal are presented. 
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1. Introduction 

The decay of a neutral kaon into two pions in a state with isospin / is described by the transition 
amplitudes 



where H w is the AS = 1 effective weak Hamiltonian and 5/ the 7T7r-scattering phase shift. In exper- 
iments it is observed that the kaon decay amplitude into two pions with total isospin / = is about 
twenty times larger than the amplitude into a state with 1 = 2, i.e. 



This enhancement, referred to as AI = 1/2 rule, remains one of the long-standing problems in 
hadron physics. Within the Standard Model, short-distance QCD and electroweak effects yield 
only a moderate enhancement. Therefore, the main contribution is expected to come from long- 
distance, i.e. non-perturbative, QCD effects or, if this is not the case, from new physics. Lattice 
QCD is the only known technique that allows to attack the problem from first principles, and pos- 
sibly to reveal the origin of the AI = 1/2 rule. 

In the low-energy regime of QCD various sources for the enhancement are possible. These include 
pionic final state interactions at around 100 MeV; physics at an intrinsic QCD scale of Aqcd ~ 250 
MeV; or physics at the scale of the charm quark, i.e. around 1.3 GeV. It remains unclear whether 
the experimental observation is the result of an accumulation of several effects, or mainly due to a 
single cause or mechanism. 

A theoretically well-defined strategy to disentangle non-perturbative QCD contributions from var- 
ious sources was proposed in Ref. [1], with the specific aim to reveal the role of the charm quark 
in the explanation of the AI = 1 /2 rule. The possibility that the enhancement is mainly due to its 
mass being decoupled from the light quark mass scale was pointed out a long time ago [2]. 

2. AI = 1/2 rule on the lattice 

In the approach of Ref. [1] the direct computation of K — > %% amplitudes on the lattice is 
bypassed by considering the K — > % and K — > vacuum transitions, which are then related to the 
physical ones by means of Chiral Perturbation Theory (ChPT) [3]. 

A crucial part of the strategy is to keep an active charm quark, such that the theory has a softly 
broken SU (4)l x SU (4)# chiral symmetry. The role of the charm quark in the dynamics underlying 
the AI = 1 /2 rule for kaon decays can then be studied by monitoring the dependence of kaon decay 
amplitudes on the charm quark mass m c . This is done in two steps: i) Set m c equal to the light 
quark masses, i.e. m u =m l i = m s = m c (GIM limit), ii) Increase the charm quark mass towards its 
physical value. 

After the Operator Product Expansion (OPE) to lowest order the AS = 1 effective weak Hamiltonian 
H w is given by 



iA ie i5 > = (J y %%)i\H w \K i 
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where V qq > are CKM-matrix elements and kf 2 the Wilson coefficients, which incorporate all the 
high-energy effects. The four-quark operators Qf are given by 

Qf = {{sy^u){uy^P.d) ± (sy M P_d)(sy M P_ B )} -(u->c). (2.2) 

Under SU(4)l the operator 2^(2]") transforms as an irreducible representation of dimension 
84(20), while both Qf are singlets under SU(4)r. In the case of a diagonal mass matrix the 
operators Qf are of the form 

<2± = {m 2 u -m 2 c ) {m d {sP+d) + m s {sP-d)} , (2.3) 

where P± = 1 /2(1 ± 75). Even though Qf do not contribute to the physical matrix elements, they 
are allowed by the underlying symmetries as a part of the effective Hamiltonian and mix with Qf 
under renormalization if m c ^m u . 

At leading order in ChPT, the ratio of amplitudes IA0/A2I is related to a ratio of low-energy con- 
stants (LECs) gf via 



Ao 
A 2 

Here, the LECs g 1 are the couplings multiplying the counterparts of the four-quark operators Qf 
in the effective Hamiltonian of the low-energy theory [3]. They can be determined by computing 
suitable correlation functions of Qf and Qf in LQCD and matching them to the corresponding 
expressions in ChPT. The matching can be performed either in the standard p-regime of ChPT, or 
in the £-regime. The advantage of the latter is that no new LECs appear at next-to-leading order, 
which a priori may allow for a better control of systematic uncertainties. 

The complicated renormalization and mixing patterns of four-fermion operators usually encoun- 
tered in lattice formulations can be avoided through the use of the Neuberger-Dirac (overlap) oper- 
ator [4], i.e. 

1 / A \ a 

D=- 1 7 = ), A = l+s-aD w , d= , \s \< 1. (2.5) 

a\ VA^aJ l + s 

Here, Dw refers to the Wilson-Dirac operator and the tunable parameter s allows to improve the 

locality properties of D. Introducing the modified quark field *P = (1 — aD/2)*P guarantees that 

the renormalization and mixing of Qf are like in the continuum theory and, in particular, that no 

mixings with lower-dimensional operators with enhanced divergences occur [5]. The combined 

use of a SU (4)-flavor symmetry and chiral fermions, thus, leaves one with logarithmic divergences 

only. 

The use of dynamical overlap fermions is computationally very expensive, and the first studies 
(e.g. [6,7]) have been carried out in the quenched approximation. Intrinsic QCD contributions to 
the enhancement can be isolated by determining gf in the theory with m u = m ( i = m s = m c . In this 
case only the figure "8"-diagram (left of Figure 1) has to be dealt with. A moderate enhancement 
is observed [7], namely 

£ «6. (2.6) 
A2 

This is not large enough to explain the experimental ratio but is already significant and cannot be 
attributed to penguin diagrams. 
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Figure 1: Diagrams to be computed: "8"-diagram (left) and "eye"-diagram (right). 

In the next step IA0/A2I has to be monitored as m c departs from the mass-degenerate limit towards 
its physical value. Thereby, the specific contribution of the charm quark to the AI = 1/2 rule can 
be investigated in detail. However, as soon as m c 7^ m u a new kind of diagram, referred to as 
"eye"-diagram (right of Figure 1), emerges and spoils the statistical signal. Its correlation function 
consists of two terms, a "color-connected" (con) and a "color-disconnected"(<i«j term which read 

C%" « (Tr {^P-S(z,x)roP-S(x,y)roP-S(y,z)^P-S u /c(z,z)}) (2.7) 
C%\ « (Tr { fyP-S(x, z) f 7oP-S(x, y) YoP-S(y, z) } Tr { y M P_S„ /c (z, z) } > . (2. 8) 

The challenge consists in the closed quark loop S u / C (z,z) with either a charm quark (c) or an up 
quark (u) running through it. By means of conventional techniques of computing quark propagators 
involving point sources, it is not possible to sample over the coordinate z and noise dominates over 
the statistical signal. More sophisticated "all-to-all" propagators have to be developed to compute 
these contractions. 



3. All-to-all propagators 

At low quark masses the numerical computation of correlation functions is hampered by huge 
statistical noise. In the spectral representation the quark propagator can be written as 

where v denotes eigenmodes and X the corresponding eigenvalues. If the quark mass m is of 
the same order as the gap, AA = — A,-, between two consecutive eigenvalues, the low-lying 
spectrum is discrete and the lowest modes have a big weight in the sum. Space-time fluctuations in 
these eigenfunctions can lead to large fluctuations in observables. Low-mode averaging (LMA) [8, 
9] has been shown to be a suitable tool to reduce these fluctuations. It amounts to separating the 
ni ow lowest eigenmodes from the rest, i.e. to truncating the sum over all eigenmodes. Consequently, 
the quark propagator decomposes into a "low"(5 ,/ ) and "high"(S*) part 

S(x,y)=S'(x,y) + S\x,y) = 1 £ V M^AM +S \ x ,y). (3.2) 

V j Ai + m 

The latter lives in the orthogonal complement of the subspace spanned by the ni ow lowest modes. 
Decomposing the quark propagator into a low-mode and high-mode propagator results in the split- 
ting of the correlation function as well. So, a common three -point function which consists of 4 
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propagators is split into 5 terms with in total 16 distinct contributions. Schematically, 

£p>pt qIIU _|_ jjlllh _|_ ^jllhh _|_ ^jlhhh _|_ ghhhh 3) 

# diagrams = 1 4 6 4 1 

where the upper indices / and h denote the number of low-mode and high-mode propagators ac- 
cording to the decomposition of eq. (3.2). The statistical signal for the correlation function is 
enhanced by exploiting translational invariance in the terms with at least one low-mode propagator. 
Since the latter is known for all space-time locations these contributions can be sampled over many 
different source points. Furthermore, to improve the statistical signal at fixed computational cost 
it turned out that instead of merely increasing n; ow it is advantageous [10] to keep n; ow reasonably 
small and construct "extended" all-to-all propagators, which allows to sample also over the high 
part of the propagator. Here, the mode itself is used as a source for an additional inversion. More 
precisely, taking the source to be the left-projected eigenmode (yoP-)v,- at the fixed timeslice t =tf, 
the solution vector S ext reads 

SrtoU= f/ = Q\s\x^t f )( Y{) P-) Vi (f,t^ (3.4) 

y 

where an average over the spatial coordinate y is automatically performed. LMA with extended 
all-to-all or conventional point-to-all propagators for the "high" part S h is crucial to obtain a signal 
in the £ -regime even for the simple case of the figure "8"-diagram. 

The "eye"-diagram, however, which is required when decoupling the charm quark mass is still 
dominated by statistical noise. To reduce its variance, stochastic volume sources (SVS) [11] and 
dilution [12] techniques are used to estimate S h stochastically. That is, low-mode averaging is com- 
bined with stochastic sources (LMA + SVS). The resulting all-to-all propagators allow to average 
over the intersection point of the 4-quark lines, i.e. the point where two quark propagators attach 
to the closed quark loop. 

To compute stochastic all-to-all propagators an ensemble of r = l,...,N T random noise vectors, 
{r]( r \xo,x)} , is generated for each gauge configuration. These source vectors are created by as- 
signing independent random numbers to all components, i.e. to all lattice sites, color and Dirac 
indices and have to obey the following two conditions 1 

{■niM) m = /i^I in {r) ) a a M = (3-5) 

1 r r=l 

^(x,x )(T7 + )^o)) src = 8 Xoya 8^8 a p8 ab . (3.6) 
Then one can invert for each of these noise vectors and obtain an estimate for the full propagator 

[K(x)(n%(y)) stc =^y), (3.7) 

where the individual solution for each of the N r noise vectors is given by 

( <j>w ):w=ii^^)( t 7 ( ' ) );w- ( 3 - 8 ) 



'Latin(Greek) letters denote color(spin) components. 
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An essential step towards reducing the intrinsic stochastic noise is the application of "dilution". 
In this work dilution is applied in spin, color and time, i.e. each of the N r noise vectors of the 
ensemble has random entries only on one spin-color component of a single timeslice with all other 
entries set to zero. 



4. Results 



A single quenched lattice is used to study the performance of LMA + SVS where volume 
sources are used for estimating S h stochastically. The bare coupling constant is /3 = 6/gg = 5.8485 
which corresponds to a lattice spacing a ~ 0. 12 fm, and a volume Va~ 4 = 16 3 x 32. The light bare 
quark mass is ami\ou = 0.02, resulting in a pion mass m n rj 320 MeV. Two charm quark masses 
are considered: am c = 0.04 = 2 x amu g ht and am c = 0.2 = 10 x amu g h t . Twenty low-modes are 
computed for each of the 120 quenched configurations. The volume sources are diluted in time, 
spin and color. Stochastic estimators are used for the loop propagator and, if required, also for the 
other propagators in order to be able to average over the position of the 4-quark line intersection. 
By applying LMA the "eye"-diagram splits into 5 distinct terms grouped by the number of low- 
mode propagators according to eq. (3.3). In the following ratios of these terms are shown. More 
precisely, the ratios are denned by 

D ( , ,x C*t{\x -zo\,\yo-zo\)±C c ™(\xo-zo\,\yo-zo\) 

Ml* " zol, bo - zol) = c 2 ( Xo )C 2 (yo) ' (4 " 1} 

where C 2 is the two-point function of the left-handed current Jo = (WjqP^). 
The overall improvement of LMA+S VS compared to LMA is illustrated in Figure 2. It reveals that 
LMA+SVS is effective for the terms consisting of 2 or 3 low-mode propagators. The variance is 
reduced significantly, most notably when the charm quark in the closed loop is heavy. In the latter 
case the absolute error of the sum of the 5 terms is roughly halved. For the term with a single low- 
mode propagator the technique of LMA+SVS shows no significant improvement; presumably in 
this case the use of multiple independent stochastic estimates for several high-mode parts increases 
the intrinsic stochastic noise and the technique deteriorates its performance. 



5. Summary and outlook 

We have reported on the progress of our ongoing project to understand the role of the charm 
quark and its associated mass scale in non-leptonic decays of kaons into two pions. When the 
charm quark mass is decoupled from the light quark masses, it is hard to obtain statistical signals for 
"eye"-diagrams. A combination of low-mode averaging and stochastic volume sources is applied 
to cure this. We observe a significant variance reduction for several contributions, even though the 
overall error remains sizable. In the next step of this project the results for the bare ratios will be 
renormalized. To this purpose the contributions of the operators Qf have to be taken into account. 
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Figure 2: Individual contributions to the ratios R± classified by their number of low-mode propagators for 
a light (top) and heavy (bottom) charm quark. Blue data points (displaced to the right for better visibility) 
result from LMA only, whereas LMA + SVS is used for the red data points (the stochastic estimate is only 
considered for a diagram if its noise is reduced, otherwise the LMA data is kept). The contribution with 
low-mode propagators, i.e. C , is not computed stochastically and C refers to the sum of all 5 terms. 
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